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Summary.  In  the  past,  a  mathematical  equation  to 
predict  the  metabolic  cost  of  standing  or  walking 
(Mw)  was  developed.  However,  this  equation  was 
limited  to  speeds  < 2.2  m  -s_l  and  overestimated 
the  metabolic  cost  of  walking  or  running  at  higher 
speeds.  The  purpose  of  this  study  was,  therefore, 
to  develop  a  mathematical  model  for  the  meta¬ 
bolic  cost  of  running  (Mr),  in  order  to  be  able  to 
predict  the  metabolic  cost  under  a  wide  range  of 
speeds,  external  loads  and  grades.  Twelve  male 
subjects  were  tested  on  a  level  treadmill  under 
different  combinations  of  speed  and  external 
load.  Speed  varied  between  2.2  to  3.2  m  •  s"  '  us¬ 
ing  0.2  ms1  intervals  and  external  loads  be¬ 
tween  0— 30  kg  with  10  kg  intervals.  Four  of  the 
subjects  were  also  tested  at  2  and  4%  incline  while 
speed  and  load  remained  constant  (2.4  m  •  \ 

20  kg).  The  model  developed  is  based  on  M»  and 
is  proportionately  linear  with  external  load  (L) 
carried  as  follows: 

Mr  =  Mw-0.5(l  — 0.01  L)(MW—  15L-850),  (watt) 

The  correlation  coefficient  between  predicted 
and  observed  values  was  0.99  ( /* < 0.0 1 )  with  SER 
of  7.7" The  accuracy  of  the  model  was  validated 
by  its  ability  to  predict  the  metabolic  cost  of  run¬ 
ning  under  different  conditions  extracted  from 
the  literature.  A  highly  significant  correlation 
( r  =  ().s> 5,  P< 0.02,  SER  =  6.5%)  was  found  be¬ 
tween  our  predicted  and  the  reported  values.  In 
conclusion,  the  new  equation  permits  accurate 
calculation  of  energy  cost  of  running  under  a 
large  range  of  speeds,  external  loads  and  in¬ 
clines. 
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Introduction 


Over  the  last  two  decades,  the  Military  Ergonom¬ 
ics  Division  of  the  U.S.  Army  Research  Institute 
of  Environmental  Medicine  has  been  establishing 
the  data  base  and  developing  a  series  of  predic¬ 
tive  equations  for  physiological  responses  of 
clothed  soldiers  performing  physical  work  in  var¬ 
ious  environmental  extremes  (Pandolf  et  al.  1986). 
Individual  predictive  equations  for  rectal  temper¬ 
ature,  heart  rate  and  sweat  loss  as  a  function  of 
the  physical  work  intensity,  environmental  condi¬ 
tions  and  particular  clothing  ensemble  have  been 
published  (Givoni  and  Goldman  1972;  Givoni 
and  Goldman  1973;  Shapiro  et  al.  1982).  These 
prediction  equations  incorporated  a  metabolic 
component  which  in  itself  can  be  predicted  by  the 
following  equation  as  originally  published  by 
Pandolf  et  al.  ( 1977): 

Mw=  l.5W  +  2.0(W+  L)(L/W): 

+  r|(W+  L)(  1.5V2 +  0.35  GV),  (watt) 

where  Mw  =  metabolic  rate  (watt),  W  =  nude  body 
weight  (kg),  L  =  clothing  and  equipment  weight 
(kg),  q  =  terrain  factor,  V  =  walking  velocity 
(m  •  s  '),  G  =  grade  (%).  This  prediction  model  of 
metabolic  cost  is  limited,  however,  to  standing  or 
walking  with  and  without  loads,  on  level  or 
graded  terrains,  but  is  not  applicable  for  walking 
speeds  above  ~2.2  m  •  s  1  or  running.  At  higher 
speeds,  the  efficiency  of  running  becomes  higher 
than  that  of  walking,  which  means  that  the  predic- 
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tion  model  of  walking  overestimates  the  actual 
energy  cost  of  running  (Ogasawara  1934;  Mar- 
garia  et  al.  1963;  Keren  et  al.  1981). 

In  the  past,  Givoni  and  Goldman  (1971)  sug¬ 
gested  that  the  metabolic  cost  of  running  (Mr) 
could  be  expressed  as  a  linear  function  of  the  pre¬ 
dicted  metabolic  cost  of  walking  (Mw)  as  follows 
(units  were  given,  originally,  in  Kcal  •  h  '): 

Mr  =  [Mw  +0.47(900-  Mw)](l  +  G/100), 

(kcal  •  h  ') 

According  to  this  model,  the  crossover  point  for 
efficiency  between  walking  and  running  is  con¬ 
stant  (900  Kcal  h  ';  1050  W).  This  prediction 
model  was  challenged  by  some  investigators  as 
overestimating  the  actual  values,  probably  due  to 
the  use  of  too  high  a  crossover  point  for  efficiency 
(Falls  and  Flumphrey  1976).  Furthermore,  the  cor¬ 
relation  coefficient  between  observed  and  pre¬ 
dicted  values  of  energy  cost  reported  by  Givoni 
and  Goldman  themselves  was  only  0.86.  This  val¬ 
ue,  while  indicating  a  good  linear  relation  be¬ 
tween  the  predicted  and  actual  values,  leaves  26% 
of  the  variance  unexplained  by  this  equation.  This 
high  level  of  unaccountable  variance  is  too  great 
to  be  acceptable  when  a  mathematical  model  is 
adopted  to  describe  this  type  of  relationship. 

The  purpose  of  this  study  was,  therefore,  to 
develop  a  mathematical  prediction  equation  for 
the  metabolic  cost  of  running.  The  intended  use 
of  this  equation  is  to  extend  the  existing  predic¬ 
tion  equation  used  for  walking,  in  order  to  be  able 
to  predict  the  metabolic  cost  under  a  wide  range 
of  running  speeds,  loads  and  grades. 


Materials  and  methods 

Subjects.  Twelve  young  fit  males  gave  their  informed  consent 
lo  participate  in  the  study.  The  subjects  had  an  average 
(mean  +  SEM)  age  of  20.6±0.6yr:  weight,  68.7±l.6kg; 
height.  172.4  +  2.1  cm:  and  body  fat,  13.1  +  1.5%.  determined 
by  underwater  weighing  according  to  Goldman  and  Buskirk 
(1461).  Their  mean  maximal  oxygen  uptake  (T„  )  deter¬ 

mined  for  uphill  treadmill  running,  measured  according  to 
the  method  of  Taylor  et  al.  (1955)  was  3.88  ±  0.09  I  •  min  ' 

( 56.8  ±  1.6  and  ml  kg  1  ■  min  '). 

Protocol  Eighteen  combinations  of  load  and  speed  (Table  I) 
were  presented  to  each  subject  in  a  randomized  order.  Before 
experiments  all  subjects  were  familiarized  with  treadmill  run¬ 
ning  with  and  without  loads.  While  running  without  loads,  the 
subjects,  wore  shorts  and  running  shoes  which  weighed 
-0.5  kg.  Loads  consisted  of  the  subjects'  clothing  and  boots 
(3.2  ±0.1  kg)  and  a  backpack  adjusted  to  give  a  total  load  of 
10,  20,  and  30  kg.  Lead  bars  were  used  for  the  external  load. 
They  were  secured  to  an  army  aluminum  backpack  frame. 
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Table  1.  Combinations  of  external  load,  speed,  and  grade 
studied 


Velocity 

Load 

Cirade 

(m  s  ') 

(kg) 

("<’) 

2.2 

0.  10,  20,  30 

0 

2.4 

0,  10,  20,  30 

0 

2.6 

0,  10,  20 

0 

2.8 

0,  10,  20 

0 

3.0 

0,  10 

0 

3.2 

0,  10 

0 

2.4 

20 

2,  4 

which  was  well  fitted  to  the  subjects.  Running  speed  ranged 
from  2.2  — 3.2  m  s-1  employing  increments  of  0.2  m  s~'. 
Each  experiment  was  10  min  in  duration  at  0%  grade  on  a 
treadmill  in  an  ambient  climate  of  22  — 25°C  and  30— 40%  re¬ 
lative  humidity.  Each  subject  completed  four  tests  per  day 
with  60—90  min  of  rest  between  tests  to  allow  for  proper  lac¬ 
tate  clearance  (Diamant  et  al.  1968).  To  determine  the  effect  of 
incline  on  the  metabolic  cost  of  running,  four  of  the  subjects 
were  tested  at  inclines  of  2  and  4%,  while  running  at  constant 
conditions  of  speed  (2.4  m  s-1)  and  load  (20  kg). 

Oxygen  uptake  (T<>,)  was  measured  by  open  circuit  spi¬ 
rometry.  The  expiratory  gases  were  analyzed  by  an  automated 
metabolic  measurement  system  (MMC  Horizon,  Sensormed- 
ic).  Measurements  were  taken  for  the  last  5  min  of  each  run  at 
15-s  intervals.  During  this  period,  the  subjects  were  at  a  steady 
state  and  T0,  did  not  vary  by  more  than  ±2.5%.  Energetic 
equivalent  of  Oj  was  presented  in  watts  (W)  based  on  the  ap¬ 
propriate  respiratory  quotient  (RQ)  during  exercise  (Consola- 
zio  et  al.  1963).  ECG  electrodes  were  affixed  to  the  subjects' 
chest  (CM5  placement)  and  heart  rate  was  determined  at  the 
same  time  as  oxygen  uptake  measurements. 

Statistics.  All  mathematical  manipulations  assumed  the  linear 
relation  between  the  metabolic  cost  of  walking  (M„)  and  that 
of  running  (Mr).  Thus  the  developing  of  the  prediction  equa¬ 
tion  was  based  on  the  transformation: 

Y  =  ax  +  b 

Multiple  linear  regression  analysis  was  used  to  compare  the 
predicted  values  M,,,.,  to  the  observations  M .  A  prediction 
equation  was  adopted  if  the  variance  of  error  ( 1  —  r- )  was  less 
than  2%  at  a  significance  level  of  p<  0.05. 


Results 

Relationship  between  M„.  and  Mr 

The  basic  transformation  to  describe  the  meta¬ 
bolic  cost  of  running  (Mr)  in  relation  to  that  of 
walking  (M„.)  is; 

Mr  =  +  (a\  +  b)  (I) 

for  V>  2.2  m  •  s  '. 

Analysis  of  all  216  observations  while  running 
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Fig.  1.  The  relationship  between  the  calculated  metabolic  rate 
of  walking  (M„)  and  the  observed  metabolic  rate  of  running 
The  dashed  line  indicates  the  mean  for  all  subjects:  the 
solid  line  is  the  line  of  identity 


or 

Mr=Mw-0.45(Mw- 1000)  (3a) 

for  F>2.2  m  •  s- 

The  effect  of  load  carriage 

Based  on  previous  observations  which  showed 
that  load  carriage  affects  the  crossover  point  be¬ 
tween  walking  and  running,  the  value  given  in  Eq. 
3  is  presumably  a  value  which  should  not  be  con¬ 
stant.  Both  the  y  intercept  (crossover  point  be¬ 
tween  walking  and  running)  and  the  slope  of  the 
regression  line  were  found  to  be  a  function  of  the 
weight  of  the  load  carried.  This  results  in  the  fol¬ 
lowing  transformation  of  Eq.  3: 

Mr  =  Mw  +  0.5(l  -0.01L)(850+15L-MW)  (4) 


on  the  level  treadmill  with  different  loads,  re¬ 
sulted  in  a  linear  relationship  between  the  ob¬ 
served  metabolic  rates  and  the  predicted  values 
for  walking  as  is  presented  in  Fig.  1.  Eq.  1  can  be 
transformed  to  be  a  function  of  Mw: 

M,  =  M.  +  (aMu+fe)  (2) 

In  this  equation,  h  is  the  crossover  point  of  effi¬ 
ciency  between  walking  and  running  and  a  is  the 
slope  of  the  line.  The  best  fit  regression  line  de¬ 
scribing  this  transformation  is: 

Mr  =  M«  +  (  -  0.45  M*  +  450)  (3) 


for  V>  2.2  m  •  s  “  '. 

The  term  (850+ 15L)  in  Eq.  4  indicates  that 
the  crossover  point  (i.e.,  y  intercept)  between 
walking  and  running  is  a  function  of  the  external 
load.  The  term  [0.5 ( 1  -0.01  L)]  indicates  that  the 
change  in  slope  of  the  regression  line  is  also  a 
function  of  the  load.  This  equation  can  be  trans¬ 
formed  to  show  that  Mw  at  a  given  speed 
( >  2.2  m  •  s  "  ')  over  estimates  the  actual  metabolic 
rate  as  follows: 

Mr  =  Mw  — 0.5(1  — 0.0IL)(Mw—  15L  — 850)  (4a)  □ 

□ 

The  correlation  coefficient  between  observed 
and  predicted  values  calculated  according  to  Eq.  — 


(wad) 


Fig.  2,  a  The  relationship  between  the  predicted  metabolic  cost  of  running  (M,,,,,)  and  the  observed  The  regression  line  for 

each  subiect  is  presented  b  The  correlation  coefficient  between  predicted  and  observed  values  of  the  metabolic  cost  for  running. 
Each  point  represents  the  mean  value  of  each  load-speed  combination  of  all  subjects 
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4  is  very  high  and  is  presented  by  the  regression 
lines  in  Figs.  2a,  2b.  Fig.  2a  presents  the  individ¬ 
ual  regression  lines  for  each  subject  and  Fig.  2b 
the  regression  line  based  on  the  mean  value  for 
each  speed  and  load  combination. 


The  effect  of  grade  on  Mr 

Eq.  4  fits  very  well  with  data  collected  while  run¬ 
ning  on  the  level.  Using  this  equation  to  predict 
the  energy  cost  of  running  on  a  graded  treadmill 
(Fig.  3),  the  mathematical  expression  was  also 
highly  accurate  without  the  need  of  a  correction 
factor.  Thus,  no  attempt  was  made  to  correct  Eq. 
4  and  adjust  it  for  grade. 

Comparison  of  predicted  values  with  results  of 
other  investigators 

Figure  4  presents  the  comparison  between  the 
predicted  energy  cost  calculated  according  to  the 
present  equation  (Eq.  4)  with  the  reported  data 
from  various  other  studies.  These  reported  data, 
which  were  extracted  from  five  different  reports 
(Margaria  et  al.  1963;  Pugh  1970;  Falls  and 
Elumphrey  1976;  Keren  et  al.  1981;  Francis  and 
Floobler  1986)  represent  experiments  over  a  wide 
range  of  running  speeds  (2.2  — 3.6  m  s'1),  grades 
(0  or  5%)  and  loads  (0  or  20  kg).  A  highly  signifi¬ 
cant  correlation  exists  between  the  predicted  and 
reported  values.  Furthermore,  all  observations, 
except  two,  fall  within  the  range  of  ±  10%  from 
the  line  of  identity  and  all  of  them  are  within  the 
range  of  ±  1 5%. 


Hr.  .1.  Metabolic  cost  of  running  on  a  2  and  4"»  graded  tread¬ 
mill  <  omparison  between  predicted  and  observed  values 


Fig.  4.  (  omparison  between  the  predicted  values  of  metabolic 
cost  for  running  and  reported  values  in  the  literature.  The 
dashed  lines  represent  the  range  of  ±  1 0°/o  from  identify  (closed 
symbols  running  on  }°a  grade I 


Discussion 

Energy  expenditure  of  walking  and  running  may 
vary  within  wide  limits  due  to  a  number  of  fac¬ 
tors,  such  as  total  weight  (body  weight  plus  exter¬ 
nal  load),  speed,  grade,  and  type  of  surface  (Gi- 
voni  and  Goldman  1971;  Pandolf  et  al.  1977; 
Keren  et  al.  1981).  In  the  past,  a  mathematical 
equation  was  established  to  predict  the  metabolic 
rate  of  walking  with  and  without  loads  (Pandolf 
et  al.  1977).  This  prediction  model  incorporated 
these  variables  and  was  found  to  be  valid  for 
speeds  slower  than  2.4  m  •  s'  '.  The  present  study 
extends  this  prediction  model  over  a  wider  range 
of  speeds  (up  to  3.2  m  -  s'  ')  with  and  without  ex¬ 
ternal  loads. 

The  equation  developed  is  similar  in  nature  to 
that  suggested  by  Givoni  and  Goldman  (1971) 
based  on  a  linear  relation  between  the  metabolic 
cost  of  running  to  that  of  walking.  However,  Gi¬ 
voni  and  Goldman  based  their  prediction  model 
for  the  energy  cost  of  running  on  a  constant  mean 
value  of  1050  W  as  the  crossover  for  efficiency  be¬ 
tween  walking  and  running.  This  value  was  ques¬ 
tioned  by  Falls  and  Humphrey  (1976)  and  did  not 
match  with  the  observations  made  in  the  present 
study.  Thereafter,  Keren  et  al.  (1981)  suggested 
that  the  crossover  point  of  efficiency  may  change 
with  load  carriage.  In  addition,  it  was  observed,  in 
the  present  study,  that  the  slope  of  the  line  is  also 
dependent  on  the  load  carried.  Thus,  in  the  pres¬ 
ent  study,  both  the  crossover  for  efficiency  and 
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the  rate  of  change  in  energy  cost,  were  trans¬ 
formed  from  constants  into  varying  factors  which 
were  dependent  on  the  weight  of  the  load  carried 
( L). 

An  important  factor  in  determining  the  energy 
cost  of  walking  and  running  is  the  placement  of 
external  load  over  the  body  (Soule  and  Goldman 
1969;  Pandolf  et  al.  1977).  Manual  transport  of  an 
external  load  rnay  cause  a  forward  shift  in  the 
body's  center  of  gravity,  which  may  alter  effi¬ 
ciency  and  in  turn  may  force  an  increase  in  energy 
expenditure.  This  is  particularly  the  case  when  ex¬ 
ternal  loads  are  carried  away  from  the  body's  cen¬ 
ter  of  gravity,  such  as  on  the  ankles  or  in  the 
hands  (Kumon  and  Belding  1971;  Robertson 
1 9S2 :  Francis  and  Hoobler  1986;  Legg  and  Ma- 
hanty  1986).  For  greatest  efficiency  and  stability, 
however,  the  load  should  be  kept  as  close  as  pos¬ 
sible  to  the  trunk  and  the  center  of  gravity  of  the 
body.  It  is  suggested  that  heavy  loads  should  not 
be  carried  using  small  muscle  groups,  but  rather 
be  carried  around  the  waist  or  in  a  backpack 
(Soule  and  Goldman  1969;  Robertson  1982;  Legg 
and  Mahanty  1985).  Therefore,  Pandolf  et  al. 

( 1977)  developed  their  prediction  model  only  for 
backpack  load  carriage,  and  the  same  attitude  was 
taken  in  the  present  study.  In  reality,  except  per¬ 
haps  in  athletic  events,  carrying  a  backpack  is  as¬ 
sociated  with  wearing  heavy  footwear  (boots). 
Thus,  the  external  load  carried  by  the  subjects 
was,  in  the  present  study,  comprised  of  their 
clothing,  boots,  and  the  backpack.  No  attempt 
was  made  to  separate  the  effect  of  the  backpack 
and  boots  on  energy  cost.  Recent  studies  imply, 
however,  that  each  100  g  increase  in  total  weight 
of  footwear  causes  a  0.7  1.0%  increase  in  energy 

cost  (Jones  et  al.  1984;  Jones  et  al.  1986;  Legg  and 
Mahanty  1986).  This  increase  in  energy  cost  due 
to  footwear  might  be  important  for  athletes  and 
especially  for  distance  runners  under  a  training 
program  when  a  higher  energy  expenditure  would 
he  more  demanding  for  the  cardiorespiratory  sys¬ 
tem;  it  will  increase  work  intensity  (  L,,  /  L(>  ) 

and  ultimately  improve  T, ,  .  Under  competi- 

ti\e  running,  lower  work  intensity,  while  running 
with  lighter  shoes,  might  be  less  fatiguing.  How¬ 
ever,  for  soldiers  who  carry  backpacks,  the  effect 
of  footwear  weight  on  total  energy  expenditure 
will  usually  be  small.  The  prediction  model  pre¬ 
sented,  being  valid  for  backpack  carriage,  might 
therefore  be  somewhat  inaccurate  when  applied 
for  predicting  the  energy  cost  of  running  with 
loads  carried  by  hand  or  around  the  ankle,  which 
recently  has  become  very  popular  in  physical 
training. 


Speed  of  running  and  grade  of  incline  are  two 
of  the  factors  which  determine  energy  cost.  Mar- 
garia  et  al.  ( 1963)  who  developed  a  nomogram  de 
novo  to  predict  the  energy  cost  of  running  use 
these  two  basic  factors  as  determinants.  Givoni 
and  Goldman  (1971)  who  suggested  an  equation 
to  predict  the  energy  cost  of  running,  based  on  the 
energy  cost  of  walking,  included  grade  as  a  factor 
also  in  the  prediction  of  M,.  The  present  predic¬ 
tion  equation,  however,  was  found  to  be  depend¬ 
ent  on  external  weight  and  was  accurate  for  pre¬ 
dicting  the  metabolic  cost  of  running,  on  the  level 
or  uphill,  without  the  need  to  correct  for  speed  or 
grade.  This  probably  derives  from  a  sufficient  re¬ 
lative  weight  already  given  to  speed  and  grade 
while  calculating  M,;  whereas  external  weight 
might  have  an  additional  effect  on  the  biome¬ 
chanics  and  energy  cost  while  running. 

The  validity  of  the  prediction  model  was 
tested  by  its  ability  to  predict  the  metabolic  cost 
of  running  under  different  conditions.  The  results 
reported  in  different  publications  (Margaria  et  al. 
1963;  Pugh  1970;  Falls  and  Humphrey  1976; 
Keren  et  al.  1981;  Francis  and  Hobbler  1986) 
were  compared  to  the  calculated  values  using  the 
basic  parameters  reported.  Bearing  in  mind  the 
difficulties  to  extract  very  accurate  data  from 
published  reports,  the  highly  significant  correla¬ 
tion  between  the  reported  and  (he  predicted  val¬ 
ues  (r  =  0.952, p< 0.02)  and  the  fact  that  almost  all 
data  points  fall  in  the  range  of  ±  10%  from  the 
line  of  identity  is  very  encouraging.  However,  a 
very  limited  data  base  exists  in  the  literature.  The 
existing  data  refers  mainly  to  tests  on  the  level  or 
at  5%  grade  with  external  loads  of  0  or  20  kgs.  The 
most  comprehensive  report,  in  this  respect,  is  that 
of  Keren  et  al.  (1981),  who  tested  their  subjects 
while  running  at  different  speeds  on  a  5%  graded 
treadmill,  without  any  load  and  with  a  20  kg 
loaded  backpack.  Their  results  are  about  5% 
higher  than  generated  by  our  prediction  model. 
These  differences,  however,  could  be  accounted 
for  bv  differences  in  technique  and  non-specific  v  ar¬ 
iability  between  laboratories.  Nevertheless,  fur¬ 
ther  information  for  a  wide  range  of  loads  and 
grades  will  be  helpful  to  further  establish  the  pre¬ 
diction  model  developed  over  a  larger  range  of 
grades. 

In  conclusion,  a  highly  accurate  mathematical 
model  was  developed  to  predict  the  metabolic 
cost  of  running.  While  other  prediction  equations 
for  running  exist,  they  are  limited  to  predicting 
energy  expenditure  for  level  treadmill  running 
only  (Van  der  Walt  and  Wvndham  1973)  or  run¬ 
ning  at  very  slow  speeds  (Shephard  1969).  The 
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present  model  together  with  the  model  developed 
by  Pandolf  et  al.  ( IV77)  enables  the  prediction  of 
metabolic  cost  for  walking  and  running  at  a  wide 
range  of  speeds,  external  loads  and  grades.  How¬ 
ever,  the  present  equation  may  be  inaccurate  for 
running  with  loads  in  the  hands  or  on  the  an¬ 
kles^ 
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